Anomalous Dimension of Dirac's Gauge-Invariant Nonlocal Order Parameter in 

Ginzburg-Landau Field Theory 
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In a Ginzburg-Landau theory with n fields, the anomalous dimension of the gauge-invariant 
nonlocal order parameter defined by the long-distance limit of Dirac's gauge-invariant two-point 
function is calculated. The result is exact for all n to first order in e = 4 — d, and for all d 6 (2, 4) 
to first order in 1/n, and coincides with the previously calculated gauge-dependent exponent in the 
Landau gauge. 



I. INTRODUCTION 

An outstanding problem in gauge theories is the con- 
struction of physical correlation functions or propagators 
of the charged matter fields. As such objects involve 
fields located at different points in spacetime, the stan- 
dard forms, expressed solely in terms of the matter fields, 
are in general not gauge invariant and, consequently, not 
physical. The principle of gauge invariance by itself does 
not yield a unique prescription, and various solutions 
have been proposed a long time ago, notably by Dirac 1 
and by Schwinger— . These proposals have been used to 
investigate important physical problems such as anoma- 
lies, quark potentials, and order parameters distinguish- 
ing the different phases of gauge theories. Recently, the 
issue has received considerable attention in the context of 
high-temperature superconductors, where massless Dirac 
fermions coupled to a dynamical gauge field were put 
forward as an effective theory for studying the unusual 
properties of the normal state of underdoped materials^. 

In this Letter we contribute to this issue by show- 
ing that different gauge-invariant proposals for corre- 
lation functions lead to different physical results. We 
do so by considering the abelian Higgs or Ginzburg- 
Landau model, which describes a great variety of physi- 
cal systems, ranging from scalar QED, superconductors 
over liquid crystals to cosmic strings and vortex lines in 
superfluidsA. The model consists of a |$| 4 -theory cou- 
pled minimally to the electromagnetic gauge field A^, 
with fi = 1, ■ ■ ■ d. Its Hamiltonian is 



H = \{d li +ieA IJ 



-m 2 \ct>t 



(1) 



where F^p = d^A v — d v A^. The scalar field cj> has 



n/2 complex components and a 0(n)-symmetric self- 
interaction with coupling constant A. The coefficients 
e and m denote electric charge and mass parameters of 



the complex <p field, respectively. The last term with pa- 
rameter a fixes a Lorentz-invariant gauge. We use mostly 
the notation of statistical field theory in d space dimen- 
sions. The results apply, however, equally to quantum 
field theory in d spacetime dimensions in the Euclidean 
formulation. 

In the following, we will work at criticality by set- 
ting m = 0. The free correlation function G{x — x') = 
(<j)(x)<ffl {x'))o of the scalar field is the Fourier transform 
ofl/fc 2 : 

d d k e kx 



G{x) 



r(d/2 



LIV 



d-2 ' 



(2) 



(27r) d k 2 Ait d l 2 

For the free correlation function D 
(A fl (x)A^(x'))o of the gauge field, we must Fourier trans- 
form 



and obtain 

D (x)-W^ 



l + cO<S„„ + (d-2)(l-a) 



(3) 



(4) 



In the presence of interactions, the expectation value 
(<p(x)<f>l (x 1 )) is an unphysical quantity since it is not in- 
variant under gauge transformations 



4>{x) 



<p(x), 



d„A{x). (5) 



In fact, it vanishes identically due to Elitzur's theorem 5 . 
A gauge-invariant correlation function was first proposed 
by Diraoi. Adapted to our purposes, it reads 

G(z - x 1 ) = (0(aO e lm S A ** J ^* )A ^ z) ^{x')) . (6) 

The average denoted by angle brackets is taken with re- 
spect to the full Hamiltonian JlJ, and the external cur- 
rent Jfj,(z) satisfies the equations 



6{z 



5{z 



d 2 J^{z) 



(?) 
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PIG. 1: Plow lines corresponding to the Dirac proposal for the 
external current J M (z) specified by Eqs. Q in two dimensions. 



where the first ensures the conservation of the external 
current in the presence of a source of strength +1 at a; 
and a sink of strength —1 at x'&. The explicit form of 



the external current (see Fig.[lJ is J^{z) 



x), where 



d d k k 



{2-K) d k 2 



V ik-z 



r(d/2 - 1) 

4 7r d/2 



(8) 

Being nonlocal, the current in Eq. JBJ is more properly 
denoted by J^(z; x, x'). At the critical point, the gauge- 
invariant correlation function J§J is expected to have the 
power behavior 



G{x) 



(9) 



with the Fisher exponent r?Gi- In the ordered phase, the 
correlation function (JgJ has the large-distance behavior 



G(x 



x<) -> \4>\ 2 , 

| SE — SC' | — »QO 



where 4>{x) is the nonlocal order parameter 

}{x) = e- ie f ddzJ '^ z - x)A ^ z) 4>{x). 



(10) 



fill 



Since J'^z) is a total derivative [see Eq. J5J], <fi{x) reduces 
after a partial integration in the Landau gauge d^A^ = 
to the local form <p(x^>. In other words, J M becomes 
invisible in this gauge and the value for r]Q\ is expected to 
coincide with the gauge-dependent result for (<p(x)<p^ (x 1 )) 
obtained in the gauge a = 0. 

The purpose of this Letter is to determine ryci to first 
order in e = 4 — d and also for all d £ (2, 4) to first order 
in 1/n. 

Note that of the two equations in @, only the source 
equation is needed for gauge invariance of This can 
also be solved by the 5-function on a line L running from 
x 1 to x. 



J^{z) = <5 M (z;L) 



dz M (s) 
as — j — 6{z 
as 



x(s)). (12) 



For a straight line L with i M (s) = x'(l — s) + sx, s £ 
[0,1], this leads to Schwinger's gauge-invariant correla- 
tion function^ 



0(z) e - ie r' d5 ^ s V(a 



(13) 



whose critical properties we studied in Ref. [l2|. The 
Schwinger construction can be thought of having all the 
external current originating from the source at x and ter- 
minating at the sink at x' squeezed into an infinitely thin 
line along the shortest path connecting the two points. 
In the disordered phase, the current lines have a finite 
line tension, and this correlation function vanishes ex- 
ponentially for large distances. Because the finite line 
tension exponentially suppresses larger loops, only a few 
small current loops are present in this phase. Upon ap- 
proaching the critical point, the line tension vanishes and 
current loops can grow without energy cost. Their pro- 
liferation signals the onset of superconductivit y 8 ' 9 ' 10 ' 11 . 

The dependence of Eq. I)13|l on the shape of L can also 
be seen more formally. We observe that a deformation 
of L is a new type of gauge transformation discussed 
extensively in Refs. 



<5 M (z;L) -> 5^{z;L' 



8^{z\L) + d v 8 ll v{z;S), 



(14) 



where <5 MJ ,(z; S) is the <5-function on the surface S swept 
out in the deformation L — » L l ' . Under this gauge trans- 
formation, the correlation function l|13|l changes by a 
nontrivial phase 



4>(x) 



-ie J d d zF^(z)S^(z;S) 



0(z). (15) 



Due to the different physical content of the correlation 
functions involved, the critical behavior to be derived 
here for Dirac's correlation function JS} will be quite dif- 
ferent from that of Schwinger's 1)1 3[l calculated in our 
previous note fl2l |. 



II. e-EXPANSION 

Perturbation theory yields via Wick's theorem, three 
perturbative corrections to (j6j to lowest order in e 2 : 



G(aj-iB') = G + T +Ti +T 2 , 



(16) 



where Tq ( i ( 2 contain zero, one, and two factors of the flow 



field J M . 



The term To is calculated by standard methods. In- 
frared divergences are avoided by evaluating Feynman di- 
agrams at a finite external momentum k. Being the only 
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scale available, k is used to render dimensionful param- 
eters such as e 2 dimensionless: e 2 = e 2 /t d ~ 4 . The result 
is the well-known gauge-dependent contributioni^ii: 



V<t> 



8tt 2 



-el 



- a 
n 



3)e. 



(17) 



The lowest-order e-expansion on the right hand is ob- 
tained by inserting for e 2 the charge e 2 = 487r 2 e/n at the 
infrared-stable fixed point, which at one loop exists only 
for n > 12(15 + 4%/l5) rj 365.9. At two loop, different 
resummation techniques suggest the existence of a fixed 
point for the physical case n = 2 15 - 16 . 

Next, we calculate the last term in Eq. I|16|l: 

T 2 (x-x') = -^G{x-x') J d d zd d z'Jp(z)D llv [z-z')Jr{z') 

(18) 

which splits in a separate scalar and gauge part. Several 
integrations by part reduce the integrals in coordinate 
spacetime to the generic form 



d d z- 



\z — x\ 

27T d / 2 



r l\d-p 



T(d/2 - l-p/2) 



p r(d/2 - l)V(d/2 - p/2) \x - x'\ d - 2 -P 



-(19) 



and we obtain, with the abbreviations = d/dx^, d 1 ^ 
d/dx' 



T 2 {x 



g4 ?r 3d/2 

d d zd d z'- 



T z (dl2-\)G(x-x') 



(20) 



1 1- 



d d zd d z'- 



2 d-4 
In the limit of small e 



4 — d, this reduces to 



T 2 (x 



e 2 — -G{x -x')M.k\x -x'\). (21) 
87H 



Comparison with an expansion of (jjjj in powers of ?7gi 
gives a contribution to the Fisher exponent proportional 
to the gauge-fixing parameter a: 



V2 



s2 



8tt 2 



(22) 



This result can be checked by considering the ra- 
tio of the correlation functions (<fi(x)<p^ (x 1 )) and 
(e _le -/" d z,/ ' j ( z ) j4 m( z )^ Adapting an argument given in 
Ref . , one can show that this ratio and, consequently, 
the combination rj^ — 772 should be independent of the 
gauge-fixing parameter a. The expressions l)17|l and l)22|l 
indeed fulfill this requirement. In our previous studj>A£, 
we found for the Schwinger correlation function as only 



difference an additional contribution to r\i independent 
of a. 

We are left with the calculation of the second, or mixed 
term in Eq. I|16|l. which reads explicitly 



TAx - x 1 ) 



; 2 / d d zd d z' 



G{x -z) — G{z - x 1 ) 



dz 



xJ iy (z')D fll/ (z - z'), 



(23) 



where the right-minus- left derivative d p = dp — dp, op- 
erates only within the square brackets. To logarithmic 
accuracy, we can writeAi 



T x [x 



eG{x 



d d z d d z' 



—G{z-x> 

dz u 



dz. 



-G{x 



J v {z')Dpv{z 



(24) 



Both terms in the square brackets give the same contri- 
bution. Proceeding in the same way as before, we find 



-2T; 



2, 



(25) 



and therefore as contribution to the Fisher exponent 

-2^72 = - ^e 2 . (26) 



47T 2 



This contribution, which is again proportional to the 
gauge-fixing parameter a, is identical to the one found for 
the Schwinger correlation function- 1 - 2 . Added together, 
we obtain for the manifestly gauge-invariant correlation 
function (J^Jl 

" A*1 = (27) 
8-n"^ n 



VGl = V<p + Vl+V2 



As expected, this result for the nonlocal Dirac order pa- 
rameter coincides with the value for t]^ obtained in the 
Landau gauge (a = 0). When a is considered a run- 
ning coupling constant of the theory (JlJ, the Landau 
gauge emerges as a fixed point a* = of the renormaliza- 
tion group 1 ^. This is a special case of the more general 
result^ that a = is always a fixed point when consid- 
ering a gauge-fixing term of the form (L^A^) 2 /2a. The 
choice L M = dp then leads to the Landau gauge, while 
the choice L M = n^, with a constant vector, leads to 
the axial gauge. 

The expression 1)2 7|l is to be contrasted with ?7gi = 
— (3/47T 2 )e 2 we 12 derived for Schwinger's gauge-invariant 
correlation function JT3J. It follows that the exponent 
1)2 7[l characterizing Dirac's correlation function is less 
negative than the one characterizing Schwinger's. The 
latter coincides with r}$ in the gauge a = — 3. For this 
value of a, the external current line connecting x to x' 
has no effect. A similar observation in the context of 
quantum chromodynamics was made in Ref. |2oT |. 
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III. LARGE-n EXPANSION 



Using the integral formula JT9J), we obtain for 772: 



The leading contribution in 1/n generated by fluctu- 
ations in the gauge field is obtained by dressing its cor- 
relation function with arbitrary many bubble insertions, 
and adding the infinite set of Feynman graphs^. The 
resulting geometric series leads to the following change 
in the pref actor of the correlation function (j3j: 



(28) 



q 2 + ne 2 [c{d)/2{d- 1)] q d - 2 ' 



where the second term in the denominator dominates the 
first for small q if d E (2,4). The constant c(d) stands 
for the 1-loop integral 



C{d): 



d d k 



(27r) d k 2 (k+p) 2 



p 2 =i 



T(2 - d/2)V 2 (d/2 
(47r) d / 2 r(d- 2) 



(29) 

where analytic regularization is used as before to con- 
trol ultraviolet divergences. To leading order in 1/n, the 
value of ?70 for d£ (2,4) read a 22 ' 2 ? 



2A-d-{d-l)[4{d-l)-da] 
n (47r) d / 2 c(d)r(d/2 + 1) ' 



(30) 



which depends on the gauge-fixing parameter a. For 
d = 4 — e, this result reduces to Eq. I)17|l obtained to 
first order in e. 

We next consider the gauge-invariant version of this. 
The term T2 in l)18|l can be evaluated as before. To ex- 
tract the dependence on ln(|a: — a;'|) it will be useful to 
replace q d ~ 2 by q d ^ 2+s in Eq. lj28]). with a dummy pa- 
rameter S, which will be taken to zero at the end. Then 
the large-n limit of the gauge-field correlation function 
becomes 



D^(q) 



2 d 



ne 2 c(d) q d - 2 +s 
or in coordinate spacetime 
8 d- 1 



, (31) 



ne 2 c{d) (47r) d / 2 r(d/2) x 



2-5 



(d - 3 + a)S^ + (1 - a 



.(32) 



Proceeding in the same way as before, we find after var- 
ious integrations by parts 



x I I d d zd d z' 



d-2 6 J 



n (47r) d / 2 c(d)r(d/2)' 



(34) 



This large-n result valid for all d £ (2,4) is once more 
proportional to the gauge-fixing parameter a, just as for 
small e in Eq. I)22|l. The result can again be easily checked 
by noting that the combination n^ — 772 is independent of 
the gauge-fixing parameter a. As for small e, the result 
for the Schwinger correlation function differs only by an 
a-independent contribution to 772- 

For the mixed term Ti, we also find for large n the 
relation l)25|l between the two contributions 7\ and T2, 
and thus 771 = —2772. This expression for 771 is identical to 
the one for the Schwinger correlation function. Adding 
the three contributions together, we arrive at 



Vgi = — 



(7 - Ad)d 



n (47r) d / 2 c(d)r(d/2 + 1)' 



(35) 



independent of the gauge-fixing parameter a. This re- 
sult, valid for all d £ (2,4), is the leading contribution 
in 1/n. As expected, it coincides with the value l)30|l for 
770 obtained in the Landau gauge (a = 0). This should 
be compared to the critical exponent for the Schwinger 
correlation function 12 l)13|l which coincides with Eq. I)3Q|I 
obtained in the traceless gauge a = 1 — d, of which 
a = —3 found in the e-expansion is a special case. In 
this gauge, where the correlation function is trace- 
less, the external current line connecting x to x 1 becomes 
invisible. Although less than for Schwinger's correlation 
function, 7701 found here is negative for small e and all 
n or for d £ (2,4) and large n. In a recent Monte Carlo 
studySi of the three-dimensional lattice model (n = 2) 
in the London limit where \(p\ = const, the large negative 
value n^ = —0.79(1) was obtained in the Landau gauge. 
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